Initializing an unmodulated spin chain to operate as a high quality quantum data-bus 
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We study the quality of state and entanglement transmission through quantum channels described 
by spin chains varying both the system parameters and the initial state of the channel. We consider 
a vast class of one-dimensional many-body models which contains some of the most relevant experi- 
mental realizations of quantum data-buses. In particular, we consider spin-1/2 XY and XXZ model 
with open boundary conditions. Our results show a significant difference between free-fermionic 
(non- interacting) systems (XY) and interacting ones (XXZ), where in the former case initialization 
can be exploited for improving the entanglement distribution, while in the latter case it also deter- 
mines the quality of state transmission. In fact, we find that in non interacting systems the exchange 
with fermions in the initial state of the chain always has a destructive effect, and we prove that it 
can be completely removed in the isotropic XX model by initializing the chain in a ferromagnetic 
state. On the other hand, in interacting systems constructive effects can arise by scattering between 
hopping fermions and a proper initialization procedure. Remarkably our results are the first example 
in which state and entanglement transmission show maxima at different points as the interactions 
and initializations of spin chain channels are varied. 



I. INTRODUCTION 

Quantum communication between different regis- 
ters/processors is a vital task in fast developing quan- 
tum technology. Using mobile particles, such as photons 
or moving electrons, for carrying information is one op- 
tion which either faces the complexity of interface equip- 
ments for different physical objects (e.g. photons and 
electrons) or needs a fine control over the bus which is 
still challenging [l], An alternative is to let the in- 
formation flow through a quantum channel, physically 
realized by a chain of permanently coupled localized par- 
ticles, exploting the dynamical properties of the channel 
itself. There are several possibilities for realizing channels 
that might serve this purpose, amongst which spin-1/2 
chains have revealed particularly suitable for transferring 
quantum information from one point to another Q-IH. 

Depending on the specific physical realization of the 
overall system, it might be easier to act on the struc- 
ture of the Hamiltonian ruling the channel dynamics or 
to prepare the channel in a specific initial state. For in- 
stance, spin chains in solid state physics represent a vast 
reservoir of possible quantum channels, characterized by 
the most diverse Hamiltonians, though with fixed pa- 
rameters 0, 0- On the other hand, initializing a spin 
chain embedded on a solid-state matrix might be a hard 
task. Quite complementary, recent progress in optical 
lattices are making a real chance out of several theoret- 
ical proposals for realizing spin chains with cold atoms 
|8l-[l3| , though with some restrictions on the structure of 
the effective spin Hamiltonians actually attainable p^ . 
Moreover, different initial states can be realized in an 
optical lattice [l^ - [l6| . and new cooling techniques [13 
also provide the possibility of reaching temperatures in 



which the magnetic phases are not disturbed by thermal 
fiuctuations and so the real magnetic ground state of the 
system becomes reachable. 

There are two essential features that characterize quan- 
tum channels made of interacting localized objects: their 
dynamics is dispersive, due to the non trivial structure 
of the many-body Hamiltonian that describes the chan- 
nel, and it depends on the initial state of the chan- 
nel itself. Dispersion is always detrimental to quantum 
information transmission, and designing a non disper- 
sive channel requires a detailed engineering of the local 
couplings [l^ [13 J which is practically hard to achieve. 
It is therefore relevant to understand up to what extent 
according to the type, parameters and the length, a ho- 
mogeneous non locally engineered spin chain is usable for 
quantum communication. In particular, and at variance 
with some recent works in which high quality transmis- 
sion is achieved independent of the system initialization 
[20I - I22I we would like to see whether or not one can im- 
prove the quality of transmission by means of a specific 
initialization. Another issue which is less studied in the 
literature (unless very few cases for the case of engineered 
chains JJ., 20] ) is the effect of Hamiltonians which do not 
conserve the number of excitations. In fact our investi- 
gation here includes a wide class of Hamiltonians which 
change the number of excitations during the time evolu- 
tion. 

We consider quantum channels realized by finite spin- 
1/2 chains with homogeneous nearest neighbor exchange 
interaction of the Heisenberg type, possibly in the pres- 
ence of a uniform external magnetic field. As for the 
initial state of the channel we consider the ferromagnetic 
state, with all the spins parallel to each other, the Neel 
state, where the spins are alternatively parallel, the state 
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FIG. 1: (Color online) Schematic picture for; 
ferring; (b) Entanglement distribution. 
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built as a series of singlets, and the ground state. In or- 
der to study the interplay between the properties of the 
channel Hamiltonian and the structure of the initial state 
in determining the quality of the transmission processes, 
we specifically deal with different Hamiltonians and dif- 
ferent initial states. We present a comprehensive study 
for the transmission quality over whole phase diagram 
of the XY and XXZ Hamiltonians with the above initial 
states. 

The structure of the paper is as follows: in Section HIl 
we introduce our scheme for quantum state transfer and 
entanglement distribution in a general language. In Sec- 
tion lllll we study free fermionic models, i.e. XY Hamilto- 
nian, while the Section IIVI is devoted to the XXZ model 
as an example of interacting systems. Finally, in Section 
|V]we comment upon our results. 



II. SET-UP 

The quantum channel consists of N spin- 1/2 particles 
sitting at sites 1 to TV of a one dimensional lattice and 
interacting through the Hamiltonian 
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where (a = x, y, z) are the exchange integrals, h is an 
external uniform magnetic field applied in the z direction, 
and (jf are the Pauli operators of the spin sitting at site /. 
We prepare the channel in some initial pure state IV'ch), 
which can be either entangled or separable with respect 
to single-spin states. An extra qubit which carries the 
information is labelled by the site index and initially 
set in some arbitrary state /Cq(0). The schematic picture 
of the system is shown in Fig. [TJa). 

At t = the interaction between the qubit and the 
channel is suddenly switched on, via 



(2) 



We use sudden switching for computational simplicity 
as the dynamics is not altered by a more realistic finite 
switching time, provided that it is small compared to the 
characteristic times set by the couplings of the Hamilto- 
nian Hi 

The total Hamiltonian H = Hch + Hi rules the dy- 
namics of the overall system, whose state at time t reads 
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where p(0) — Po(0) ® |V'ch)('0ch| and h = \. The den- 
sity matrix of the qubit N , namely p„ (i), is obtained by 
tracing out the qubits 0, 1, ...,N — 1 from p{t). On the 
other hand, one can define the superoperator £t which 
maps the initial density matrix of the qubit 0, i.e. p„{0), 
to the density matrix of qubit N at time t, i.e. Pf^it), 
according to 



p,[t) = £t[pM\- 



(4) 



When the transmission of a generic pure quantum state 
I'i/'o) is in order, according to the scheme depicted in 
Fig. [TJa), we set Po(0) = |V'o)(^ol ^^^^ quantify the qual- 
ity of transmission by the fidelity {tp^ |p„ (i)IV'o) between 
the initial state |'0o) and the final state fact, 
the quality of the channel is better assessed by the aver- 
age fidelity / (i'0o(''/'olPN(^)IV'o)' where integration is over 
the surface of the Bloch sphere and dil^^ the correspond- 
ing measure. There are essentially two mechanisms that 
cause average fidelity to deteriorate during a transmission 
process of the type we are describing: dispersion, which 
is a collective phenomenon due to the channel being an 
interacting many-body system, and local rotations. As a 
matter of fact, the average fidelity does not distinguish 
between bad transmission, i.e. dispersion of the state all 
over the chain, and good transmission with an extra ro- 
tation during the dynamics; on the other hand, the latter 
can be safely handled with an extra unitary operation on 
the site TV, thus leaving dispersion as the only destructive 
effect in transmitting quantum states. 

Assume a unitary operator R, that does not depend on 
the state IV'o), is found such that the average fidelity for 
the rotated final state R^Pj^{t)R equals the maximum 
attainable value through the specific channel, then one 
could get a quantitative estimate of the dispersiveness of 
the transmission, which might be a very useful tool for 
characterizing the channel suitability for state transfer 
processes. Aiming at such goal, we introduce the Optimal 
Average Fidelity (OAF) 

F{t) = max / dyjMR^St Uo){A\] , (5) 

where maximization over R guarantees that the effect of 
local rotations is removed. It is of absolute relevance 
that, as shown in Appendix A, the OAF can be deter- 
mined explicitly in terms of the superoperator £t 
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+ j^{mi+m2+ sign(det(7\/J^)) TO3) (6) 
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where mi(t) are the singular values, in decreasing order, 
of the matrix 



(7) 



where (j„, n = 1,2,3 are the Pauli matrices. The Opti- 
mal Average Fidelity, once computed via Eq. ([5]), gives 
a quantitative indication about how well a channel be- 
haves, as far as the pure states transmission is concerned. 
Notice that the above expression has the very same form 
of the maximal fidelity with respect to maximally entan- 
gled states [131 , which, on the other hand is a completely 
different measure (an entanglement measure, in fact). 

When entanglement distribution comes into play, 
mixed state transmission must be considered, and other 
strategies are necessary. Let us prepare qubit in a 
maximally entangled state with an isolated qubit 0', see 
Fig-lUb). The dynamical evolution cause the mixed state 
of qubit to be transmitted to qubit N , thus generating 
entanglement between qubit 0' and N. We can quantify 
the quality of such transmission by the amount of entan- 
glement between qubits 0' and N at time t. For a generic 
spin chain, when qubits and 0' are initially in the max- 
imally entangled state |$+) = (|00) |11)) /V2 one can 
write 



P,,Jt) = {I^£,) [|$+)($+|] 



(8) 



where, / represents the identity map. In the following, 
we use concurrence C(p) as an entanglement measure 
[25j for quantifying the amount of entanglement shared 
between qubits 0' and N. Notice that different choices 
of initial maximally entangled state would give the same 
result, as the set of maximally entangled state can be 
obtained from |$^) through a local unitary operation in 
0', which is isolated. 

Consistently with the transfer process, our scheme im- 
plies the existence of an arrival time when both the op- 
timal average fidelity F{t) and the concurrence C{t) = 
C{p^,^[t)) get their maximum value. In fact, due to the 
finite size of the chain, the information travels from to 
N and viceversa multiple times, and the above quantities 
displays multiple peaks during the dynamics. Through- 
out this paper we concentrate on the first peak, whose 
position defines the arrival time t = t*, as in a practical 
situation waiting for longer times is unwise due to the 
effect of decoherence. 



into a free fermionic system, which makes valuable an- 
alytical results available. Despite being mapped into a 
non-interacting system, the model in the infinite N limit 
has a rich phase diagram featuring a quantum phase tran- 
sition [26| ai h = 1 and, as far as the entanglement prop- 
erties are concerned, the divergence of the entanglement 



range when approaching the curve h 
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1, where 



pairwise entanglement vanishes [27.- - .-29.] . Moreover, its 
peculiar non equilibrium dynamics has been studied in 
the framework of dynamical entanglement sharing [30l |. 
with periodic boundary conditions assumed. 

For diagonalizing this Hamiltonian one first maps 
spin operators tr^ = (cr^ ± iCT^)/2 to fermionic oper- 
ators through the Jordan- Wigner transformation ci = 



n 



1 

n=0 ' 



-(T^)(T, where {q,c;'} = and {cuc],} 
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as can be easily proven. The fermionic Hamiltonian can 
then be diagonalized with the procedure described in [3l| , 
of which we give a short summary in appendix B. The 
resulting diagonal Hamiltonian is 
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(10) 



where the diagonal fermionic operators are obtained via 
the Bogolubov transformation 
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In the here considered case of finite N and open bound- 
ary conditions the analytical expressions of the energies 
Ek and of the matrices P, Q for finite 7 and h are com- 
plicated [3^ . but they can be determined numerically as 
explained in appendix B. 

The channel can be conveniently characterized in terms 
of its dynamical behaviour, which is investigated in 
the Heisenberg representation. The fermionic operators 
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Uit) =F^e^**^ P - 
W{t) =P'^e-'*^ Q + Q^ 
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III. FREE FERMIONIC SYSTEMS: XY 
HAMILTONIAN 

In this section we consider the XY model defined by 
Eq. (dJ with 

1 -I- 1 — '7 

J. = J^^, Jy^J^^, Jz^O, (9) 

where, J is the exchange coupling and 7 is the anisotropy 
parameter. This model is exactly solvable, as it turns 



E is the diagonal energy matrix with elements Ek, and 
the inverse transformation Eq. (|B7p together with the 
identity e+^^'d^e"*^* — e~*'^'=(ife have been used. From 
the dynamical evolution of c/(t) that of the last-spin op- 
erators follow 
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where conservation of the parity n/^o(~'^f ) used 
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We prepare the qubit in the density matrix which is 
parametrized by the vector n(0) of the Bloch sphere via 
Po(0) = + ^^(0) • i?)/2. Correspondingly, the time evo- 
lution of the state of qubit N, pN(t) = {I + n'{t) ■ (t)/2, is 
parametrized by the vector n'„{t) = {'^{0)\a%{t)\^{0)) , 
which is determined from ([T5| . Assuming that the chan- 
nel is initialized in some state with constant parity p, 

i-e- rifcLi i-'^k) IV'ch) = (-l)^|'0ch), our quantum chan- 
nel (HI is described by the following afRne transformation 
which maps vector n(0) of the Bloch sphere of qubit to 
vector n'{t) of the Bloch sphere of qubit N: 
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where u{t) = pmit)], w{t) = \Wmit)\, K{t) = 
aigpNoit)] + {v + 1)71", 4>w{t) = - aig\WNo{t)] + pn, the 
rotation matrix is 
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and A{t) is the real function 
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Eqs. p6|) show that the evolution in the xy-plane in- 



volves the two rotations 5*^ 



SI 



as well as the 



shrinking towards the center of the Bloch sphere embod- 
ied in u{t) — w{t) and u{t) + w{t). Notice that none of 
the quantities involved depends on the initial state |V-'ch), 
thus relating the dynamics in the xy-plane only with the 
phase-diagram j-h of the model. In fact, the only depen- 
dence on the initial state is in the quantity A{t) which 
uniquely affects the shift in the z direction, and repre- 
sents the interference of \ipo) with IV'ch) during the evo- 
lution. Notice that in the XX case, where Wij{t) = 
(see appendix C), the dynamics (fTB)) realizes a general- 
ized amplitude damping channel 33]. 

The OAF and the Concurrence can be calculated once 
the superoperator £t is known. The elements of the su- 
peroperator £t can be read from 



k,l=0,l 



From the elements {i\£ {\k) {l\)\j) we can construct the 
Choi matrix Cs (t) 



mc£{t)\ij) = m{\k)mj) 



(20) 



which completely characterizes [3J| both the superoper- 
ator £t and the state (|8]). By explicit calculation we have 
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(21) 



From the Choi matrix, using the procedure described 
in appendix A, we can compute the OAF: the singular 
values of Meit) are 2\u{t) + w{t)\, 2\u{t) — w{t)\, and 
2\u{tf - ■w{tY\, from which 



1 1 



1 



F^^ {t)^- + - \u{tY - iD{tY\ + - max{u(t), u.(t)}. 

(22) 

This is a remarkable result as it is fully independent of the 
initial state of the channel (i.e. the parameter A{t)) and 
depends only on the Hamiltonian parameters. Moreover, 
the rotation that maximizes the average fidelity is found 



to be 



R = 



for u > w 



(23) 



In the XX case the effect of the magnetic field on the dy- 
namics is only in the phase 0^, as shown in appendix C, 
and therefore one can always choose the magnetic field 
such that at the arrival time t*, i.e. when fidelity peaks, 
4>u{t*) = 0. On the other hand, in the XY case the dy- 
namical quantities depends on /i in a complicated way, 
making the explicit rotation R, given in Eq. (j23p . neces- 
sary. 

From Eq. we see that the larger the difference be- 
tween u{t*) and wit*) the larger the OAF. In particular. 
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FIG. 2: (Color online) (a) vs. 7 and h; the peak is 

located at 7 = 0.7 and h = 1. (b) w(t*) vs. 7 and /i. Both 
figures are for A'^ = 50. 



as seen in Fig. [21 wc find that, whenever the fidehty is 
large, we see u{t*) ^ w{t*), so the qualitative behavior 
of the OAF is the same of u{t*). Moreover, in the re- 
gion < 7 < 1 and 1 < < 2, u{t*) is large, taking 
its maximum value for 7 = (XX case) and for 7 = 0.7 
and h = 1. This means that in the 7-/1 phase-diagram 
the line 7 = and the point (0.7, 1) set the best possi- 
ble Hamiltonian parameters, corresponding to the least 
dispersive channel. 

The scaling of {t*) for increasing length N is 
shown in Fig. [3] where it is clear that for the best pa- 
rameters (7 = 0, /i = and 7 = 0.7, h = 1) the OAF 
decreases very slowly, and it is greater than the classical 
value (2/3) even for chains up to = 240. This can be 
proved in the XX case thanks to the analytical results of 
appendix C. In fact, using Eq. (jC3|) . we can show that 
the solution of equation F^^ [t*) = 2/3 is A^ = 240, in 
excellent agreement with Fig. [31 The weak dependence 
of F^^{t*) on A^ for the best parameters strengthens the 
statement that these indeed define the least dispersive 
channel, no matter the length. Conversely, for non op- 
timal parameters (for example 7 = 1, /i = 1) the OAF 
decreases quickly and becomes lower than the classical 
threshold value 2/3 for chains longer than A^ — 32. 

Let us now consider the entanglement distribution. We 
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FIG. 3: (Color online) Scaling of the OAF at the arrival time 
t* versus length N, for different parameters of the Hamilto- 
nian. 



prepare the qubit pair O'-O in a pure entangled state, 
and let qubit to be connected with the chain, while 
keeping 0' isolated, so as to gen erate an entangled pair 
O'-A^. We use the concurrence [2^ to quantify the amount 
of entanglement in Pg,j^{t). 

The density matrix of the qubit pair O'-A^, Eq. ([8]), is 
given by the Choi matrix as p^,^ = ^Cf. The concurrence 
of this state is 

C(t)=max{0, C{u{t),w(t)), C{w(t),u{t))} (24) 

where, 

Cix, y) = x- v/(y2 + A{t)){l - .t2 - A{t)). (25) 

One can see that, at variance with the OAF, entangle- 
ment distribution depends on the initial state as it is a 
function of A{t)] this is due to the fact that during the 
dynamics the initial state of qubits and 0' interferes 
with the initial state of the chain and deteriorates the 
quality of the transmission. 

In the following, we investigate different initial states 
of the channel in order to find out to which state a bet- 
ter quality transmission might correspond. In particu- 
lar, we will refer to two fully separable states, namely 
the ferromagnetic state, with all the spins aligned along 
the z direction, e.g. |0, 0,...,0), and the Neel state, 
with neighbouring spins antiparallel to each other, e.g. 
|0, 1, 0, 1, 0, 1). Moreover, we will also study two dif- 
ferent entangled initial states, namely that defined by a 
series of singlet states, and the ground state of the chan- 
nel Hamiltonian. 

Let us first consider the XX (7 = 0) model, so as to 
exploit the analytical expressions available (see appendix 
C). We can prove that the concurrence achieves its max- 
imum value, i.e. u(t), when IV'ch) is initialized in a ferro- 
magnetic state. In fact in this case, since Wij{t) = 0, it 
is 

N 

A{t) = J2 UNj{trUm{t){^Pch\c]ci\^P,i,). (26) 
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FIG. 4: (Color online) Entanglement versus 7 in a chain of 
length A'^ = 20 for different initial states while the magnetic 
field takes: (a) h — and (b) h = 0.5J. 



which is equal to (1) when ji/'ch) consists in a tensor 
product of down (up) spins. For other initial states, in 
which < A{t) < 1, the concurrence is lower than u{t). 

In Fig. HJa) we plot the concurrence as a function 
of anisotropy 7 when h = for different initial states. 
As the figure clearly shows increasing the anisotropy de- 
creases the quality of transmission. 

In the limit of 7 — > 1 the Hamiltonian Hch becomes 
Ising-like which has a poor transmitting quality. As one 
can see, in Fig. |H[a), in the absence of magnetic field 
the ferromagnetic initial state always gives the highest 
entanglement; in particular, when anisotropy 7 is small 
the difference between this initialization and the others 
is evident. 

One may improve the poor ability of entanglement dis- 
tribution in highly anisotropic chains (large 7) by switch- 
ing on the magnetic field. This is clearly seen in Fig.Ul^b) 
where we set h — 0.5. Furthermore, we notice that the 
field does not essentially affect the transmission for small 
7 and, quite surprisingly, makes the ground state the best 
possible initial state for strongly anisotropic chains. 

In order to better understand the role of magnetic field 
we consider the transmitted entanglement for different 
initial states as a function of h. In Fig. [5] we plot the 
concurrence for 7 = 0.7 and 7 = 1 and see that there 
exists a value of the field, slightly depending on the initial 
state of the chain, above which the transmission becomes 
possible, even for large anisotropics. 

The existence of an exact solution for the XY Hamil- 
tonian allows us to study the quality of transmission for 
very long chains. In Fig. [5] we plot C(t*) as a function of 
N for different initial states, which evidently differentiate 
the entanglement transmission through long chains. In 
particular, in Fig. [Sla), we see that in the XX chain the 
ferromagnetic initial state not only gives the highest con- 
currence amongst the different initializations but it also 
provides the best scaling with A^. In Fig. EUb) we plot 
C{t*) as a function of A^ for 7 = 0.7 and h = 1, i.e. for 
the parameteres that defines the less dispersive XY-like 




FIG. 5: (Color online) Entanglement versus /i in a chain of 
length A'' = 20 for different initial states while the anisotropy 
parameter takes: (a) 7 = 0.7 and (b) 7=1. 
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FIG. 6: (Color online) Scaling of the obtained entangle- 
ment at the arrival time versus length A'^ for different initial 
states: (a) isotropic XX Hamiltonian (7 = and h — 0); (b) 
anisotropic XY Hamiltonian (7 = 0.7 and h = 1). 



channel (see Fig. (2). At variance with the state trans- 
mission case, where for such parameters the transmission 
quality is as high as in the XX case (see Fig. [3]), during 
entanglement distribution through an XY chain we can 
not avoid the interference A(t) by properly choosing the 
initial state of the chain, and a strong dependence on the 
length A^ appears; Fig.lSJb) in fact shows that this gives 
rise to a significant lowering of the transferred concur- 
rence. 

Results for the Ncel state and the series of singlets 
are found to be very close to each other (therefore only 
those for the former state are plotted): This shows that 
the inherent entanglement in the initial state has a very 
little effect when a state is attached at one end of a spin 
chain. 
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FIG. 7: (Color online) Transferring properties of a XXZ chain 
with the length of A'' = 20 in its whole phase diagram for 
different initial states: (a) optimized average fidelity at the 
arrival time t* (the black line indicates the classically obtained 
fidelity 2/3); (b) obtainable entanglement at the arrival time 
t*. 



IV. INTERACTING SYSTEMS: XXZ 
HAMILTONIAN 

After studying the effect of the initial state on the 
transmission quality of the channel in free fermionic sys- 
tems, interacting models must also be considered, not 
only because they represent the large majority of many- 
body systems, but also because they are usually charac- 
terized by an extremely rich, though often difficult to be 
physically deciphered, phenomenology. 

Amongst the interacting Hamiltonians the XXZ spin 
model here discussed, and defined by Eq. ([Ij with 

*=.. = |, J,^f. (27) 

is known to describe many real systems and compounds, 
thus playing an essential role, especially in one dimen- 
sional physics 0, [ill • 

This model has a very rich phase diagram: A < — 1 is 
the ferromagnetic phase with a simple separable ground 
state with all spins aligned to the same direction. For 
— 1 < A < 1 the Hamiltonian is gapless and the region is 
called XY phase. A > 1 defines the Ned phase, where 
the spectrum is gapped and a finite staggered magneti- 
zation arises. In the Ising limit A ^ 1 the ground state 
is the Neel state. 

Dynamical transmission in interacting fermionic mod- 
els is very different from the free fermionic case, as it 
results from a complex combination of many different ef- 
fects amongst which the scattering between interacting 



excitations and the existence of localized states. As a 
matter of fact, there is no transmission process which is 
generally independent on the initial state of the channel. 
In particular we find that, depending on the value of A, 
the best transmission processes correspond to different 
initializations of the chain. 

Let us first discuss the state transfer process. In 
Fig. [TJa) we see that even the OAF strongly depends on 
iV'ch): for A > the ground state clearly gives the best 
possible initialization, while for A < the phenomenol- 
ogy is much more complicated. When considering the 
entanglement distribution in Fig. the situation is 

somehow reversed, with a specific state, namely the fer- 
romagnetic one, granting the best possible initialization 
for A < 0, and a more complex scenario for A > 0. 

The general phase diagram in Fig.[7]embodies the com- 
plex balance between the effect of interference, which can 
be varied by acting on the initial state of the chain, and 
the dispersiveness of the channel which depends on A. As 
we already noted, in the XXZ case also the OAF is af- 
fected by interference, and the overlap of F{t*) observed 
in Fig. [Zlja) in the A = point is due to its correspond- 
ing to a non-interacting model, for which, as proved in 
section Hm the OAF does not depend on the initial state. 
On the other hand, there is no A value where the entan- 
glement distribution is independent on the destructive ef- 
fects of interference. Therefore, the only mechanism for 
removing such effect is that of choosing a ferromagnetic 
initial state, whose dynamics can be resolved in the sin- 
gle particle sector [3] where interference does not occur. 
The absence of interference is the reason why for A = 
the C{t*) obtained with the ferromagnetic initial state is 
by far larger then the others. Finally, in Figs. [Tfa) and 
(b) we see that, given the initial state, F{t*) — 2/3 and 
C{t*) have a similar behaviour as a function of A, with 
a shift in C{t*) for the FM initial state, consistent with 
the exact analysis given above for the XX case. 

Indications about the effects of dispersiveness of the 
channel on the transmission processes can also be ex- 
tracted: the ground state is seen to minimize such effects 
for whatever A, though this implies a better entangle- 
ment distribution only for A > 0.5 where interference 
probably plays a minor role leaving dispersion as the 
main destructive effect. 



V. CONCLUSION 

In this paper we have studied the quality of state 
and entanglement transmission through quantum chan- 
nels described by spin chains varying both the Hamilto- 
nian parameters and the initial state of the channel. We 
have considered a vast class of Hamiltonians, including 
interacting and non-interacting fermionic systems, which 
contains some of the most relevant experimental realiza- 
tions of one dimensional many-body systems, both in the 
framework of solid state physics and in the realm of cold 
atoms in optical lattices. 
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We find that if a free-fermionic model is available and 
an XY-like spin Hamiltonian can be effectively realized, 
then the best possible tuning of the parameters is that 
corresponding to the XX model with a ferromagnetic 
initial state, both for state and entanglement transfer, 
whose quality stays surprisingly high even for chains as 
long as iV ~ 240. In the anisotropic case, a state transfer 
of the same quality as that attained at 7 = 0, is ob- 
tained for 7 = 0.7 and h — 1: referring to the framework 
developed in Ref. [221, i^^^' that the relevant exci- 
tations lie in the linear zone of the dispersion relation 
and the resulting dynamics is essentially dispersionless. 
Moreover, good results for both state and entanglement 
transmission are found in a wide range of the parame- 
ters 7 and h, providing the channel is initialized in its 
ground state. When an interacting XXZ model with a 
specific A is at hand, one has to choose whether to op- 
timize the state or the entanglement transmission, since 
these goals are obtained with different initial states. In 
fact, we find that the optimal average fidelity is more 
sensitive to the dispersiveness of the channel, while the 
entanglement distribution is more sensitive to the inter- 
ference with the initial state of the chain. As a matter 
of fact the former gets its maximum in the antiferromag- 
netic isotropic (A = 1) channel initialized in its ground 
state, while the latter is maximized by an XX (A = 0) 
channel initialized in a ferromagnetic state. 

Our analysis show that the fidelity and the entangle- 
ment do not necessarily quantify the quality of quantum 
communication in the same way. Namely, highest en- 
tanglement transfer can occur along a spin chain which 
is different from that giving the highest average fidelity. 
In fact, to the best of our knowledge, these results are 
the first example in which state and entanglement trans- 
mission show different features due to the different role 
played in such processes by dispersion, essentially set 
by the parameters of the Hamiltonian, and interference, 
which explicitly depends on the initial state of the chan- 
nel. However, when we have higher entanglement one 
can always purify/distill entanglement using local oper- 
ations [35| and subsequently use it for teleportation and 
eventually end up with a higher fidelity. 
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Appendix A: Optimal average fidelity 

In this section we derive an explicit expression for cal- 
culating the optimal average fidelity (OAV) defined in 



Eq. dSJ. Using the results of Ref. ^ 

max ($+!(/ (gji?!'^) n$+)($+|l (/®i?)|$+> ^ 

Reu{2) 

max (V-meK/^f) [|*+)($+|] l^mc), (Al) 

since for varying unitary matrix R, generates 

every maximally entangled state IV'mo)- Moreover, in [23| 
it was proved that 

max(V'mo|p|'0mo) = 7 (1 + ^1 + ^2 - sign(det(T))i3) 

(A2) 

where T„in — Tr [pcTm <S5(Tn\ and ti are the singular values 
of T, i.e. the square root of the eigenvalues of T^T, in 
decreasing order. The explicit expression of equations ^ 
and ^ follows directly from Eqs. (PT]) . (|X2|) and for the 
properties of Choi matrix = {I®£) [2|<i>+) (<f>+ 1] which 
can be easily proved from its definition ((20|) . In fact, 

Ml,^ = Tr[a„£((7„)] = Tr [C^gI ® g^] = 

= (-!)"+! Tr[C^a„ $5 a„] . (A3) 

Appendix B: Diagonalization of the XY chain 

The total XY Hamiltonian is 

N-l N 

^ jY,{a+G^^^+-ia+a+^^+h.c.} + hY,al (Bl) 

fc=0 fc=0 

which, by substituting the spin operators with their 
fermionic counterparts takes the form 

H=Y. {A^kici + ^{clBkic\ - CkBkiCi)}, (B2) 

where, A is a symmetric matrix with elements A^i = 
J{Sk,i+i + 5k^i-i) + h5k,i and B is an anti-symmetric ma- 
trix with elements Bki = j{Sk,i-i — Sk,i+i)- The above 
quadratic Hamiltonian can be diagonalized in the form 
([TU]) using a Bogoliubov transformation (ITTT) , and the di- 
agonalization process reduces to finding the energies Ek 
and the matrices P and Q. From Eq. (flUl) 

[H,dk]^Ekdl. (B3) 

and then, by using Eq. (jB2l) and ((TT|) the conditions read 

AP^ + BQ'^ = +P'^E 

AQ^ + BP'^ = -Q^E (B4) 

where, i? is a diagonal matrix with elements Ej^i — 
Ek5k,i- To ensure that the transformation pT|) be canon- 
ical and invertible the matrices have to satisfy these con- 
ditions 

PP^ + QQ^ = P^P + Q^Q^In+1 
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PQ^ + QP^ 



P^Q + Q^P = 0, 



(B5) 



which can be simplified by defining the new matrices a = 
P + Q and (3 = P -Q. In fact, Eqs. dBSj force a and /3 
to be orthogonal matrices. Moreover, since {A + B)"^ = 
A — B, Eqs. (jB4[) simplifies into a single equation 



A- B ^a^ EI3 



(B6) 



which is the singular value decomposition of the matrix 
A—B. The desired matrices P = {a+(3)/2 and Q = (a — 
/3)/2 are computed accordingly and the diagonalization 
process is completed. 

For completeness we show also the inverse transforma- 
tion that easily comes from (|B5p 



where /3 = 2Jt. In the above equation we used 
the Jacobi- Anger expansion and some properties of the 
Bessel function J'n [331 • The approximation consists in 
neglecting the Bessel functions of order in{N + 2), with 
m > 2, since they contribute only after times of order 
In fact, one can show that at the transmission 
*, i.e. the time when Uj^g{t) takes its first peak, 
l3* = 2Jt* ~ N-^{N/2)k Using the properties of Bessel 
function [331 



niN 
2J ■ 

time t 



N 



ck = ^PiA + QiA^ 



1=0 
N 



4 = Y^pA + q^a- 



i=0 



(B7) 



u{t*) 



2i_ 



Ai(0 



5N 



(3eAi'(0 + 22Ai(0) (C2) 



Appendix C: Analytical evaluation of u{t) in the XX 
model 

In the XX model we have Ek = 2Jcos(^^|^) + 2h, 



Pk7i ~ 



sin C''-''"^^^^^^-' ) and Qkn = 0, making 

Wknit) = 0. Since the magnetic field causes only a con- 
stant shift in the dispersion relation, in the following we 
set h = Q. In this case 



where Ai(^) is the Airy function. It can be proved [J| 
that the maximum of u{t*) is reached for ^ — —1.019 
and thus 



N 



PkoPkN 



J2 



m—1 

4 

W 



N 

E< 

fc=0 



■Km{k -|- 1) 
N + 2 

{JnW + 2JN+2W) + JN+iifi)) 



PkoPkN 



N + 2 



JN+2W) 



(CI) 



2.700 4.804 



iV3 



N 



(C3) 
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